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We study the renormalization flow of axion electrodynamics, concentrating on the non- 
perturbative running of the axion-photon coupling and the mass of the axion (like) particle. Due to 
a non-renormalization property of the axion-photon vertex, the renormalization flow is controlled 
by photon and axion anomalous dimensions. As a consequence, momentum-independent axion 
self-interactions are not induced by photon fluctuations. The non-perturbative flow towards the 
ultraviolet exhibits a Landau-pole-type behavior, implying that the system has a scale of maximum 
UV extension and that the renormalized axion-photon coupling in the deep infrared is bounded from 
above. Even though gauge invariance guarantees that photon fluctuations do not decouple in the 
infrared, the renormalized couplings remain finite even in the deep infrared and even for massless 
axions. Within our truncation, we also observe the existence of an exceptional RG trajectory, which 
is extendable to arbitrarily high scales, without being governed by a UV fixed point. 



I. INTRODUCTION 

The existence of a fundamental pseudoscalar particle is 
strongly motivated by the Peccei-Quinn solution to the 
strong CP problem pQ. The axion, being the pseudo- 
Nambu-Goldstone boson of a spontaneously as well as 
anomalously broken axial U(l) symmetry [2], receives 
a small coupling to electromagnetism in generic mod- 
els [31 H], very similar to the neutral pion in QCD. As 
a consequence, the effective theory below the QCD scale 
contains photons and axions as light fundamental degrees 
of freedom. Their interaction is governed by a dimension 
five operator with a coupling being inversely proportional 
to the (presumably high) scale of Peccei-Quinn symmetry 
breaking. For reviews, see, e.g., [5]. 

The resulting effective theory, axion electrodynamics, 
actually has a wide range of applications. For instance, it 
also occurs in the context of macroscopic media exhibit- 
ing the magneto-electric effect [6] , as well as in any hypo- 
thetical theory beyond the standard model with axion- 
like degrees of freedom [7[ [8], and, of course, for the de- 
scription of pion decay into two photons. 

In the present work, we take a more general viewpoint 
on axion electrodynamics and study its renormalization 
properties in the effective framework with pure photon 
and axion degrees of freedom. There are several motiva- 
tions for this study: first, experiments actively searching 
for the axion cover a wide range of scales: solar observa- 
tions such as CAST [3] or the Tokio helioscope jTU] search 
for axions emitted by the sun through the Primakoff 
process with a typical keV momentum scale Even 
higher momentum scales are probed by the stellar cool- 
ing rate in the Helium burning phase of horizontal branch 
stars and other astrophysical considerations as reviewed 
in p~2l414j . with hints of possible observable effects even 
at the TeV scale [TS]. By contrast, purely laboratory 
based experiments as the classical light-shining-through- 
walls (LSW) [P6TI22] . see also [23] for an overview, or 
polarimetry set-ups pTJ EH [25] using optical lasers work 
with momentum transfers on the /zeV scale. Searches 



which are sensitive to axion (like) dark matter operate 
at similar scales [26] . A sizable running of the couplings 
over these nine orders of magnitude would have severe im- 
plications for the comparison of experimental results |27j . 
Even though the coupling is expected to be weak, pho- 
ton and axion fluctuations could potentially lead to siz- 
able renormalization effects, because of their small mass. 
In particular, photon fluctuations strictly speaking never 
decouple as their masslessness is granted by gauge invari- 
ance. 

Another motivation is of more conceptual type: even 
though QCD-type models of axion electrodynamics do 
have a physical cut-off, approximately given by the scale 
of chiral symmetry breaking ~ 0(1) GeV, it remains an 
interesting question as to whether axion electrodynam- 
ics could be a self-consistent fundamental quantum field 
theory. Within a nonperturbative context, the negative 
mass dimension of the axion-photon coupling does not 
prohibit the existence of the theory as a fundamental 
theory, it only precludes the possibility of the coupling 
becoming asymptotically free. If the theory can only 
exist as an effective theory - as naively expected - its 
breakdown in the ultraviolet can put restrictions on the 
physically admissible values of the coupling and the ax- 
ion mass, much in the same way as the Higgs boson mass 
in the standard model is bounded from above by UV 
renormalization arguments, see, e.g., [28] . 

Since the axion-photon-interaction as a dimension five 
operator is perturbatively non-renormalizable, a nonper- 
turbative treatment is required. This work is thus based 
on the functional renormalization group (RG) allowing to 
explore a controlled approach to the potentially strongly 
coupled UV. 

The motivating questions are answered by the present 
study with little surprise, though with interesting lessons 
to be learned from field theory: in the physically relevant 
weak-coupling parameter regime, the renormalization ef- 
fects remain tiny despite the presence of massless fluctu- 
ations. Furthermore, no indications for generic UV com- 
pleteness are found beyond perturbation theory. Never- 
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theless, the present study yields an interesting example 
of renormalization theory revealing several unusual prop- 
erties: we identify a nonrenormalization property of the 
axion self-interactions and show that the RG flow to low- 
est order is determined solely by the anomalous dimen- 
sions of the photon and the axion field. In addition to 
explicit solutions to the RG equations, we also identify an 
exceptional RG trajectory which exemplifies a new class 
of potentially UV-controlled flows. 

Despite the non-decoupling of massless modes in the 
deep infrared, our RG flow predicts that the physi- 
cal couplings reach finite IR values. Explicit solutions 
demonstrate, how the massless modes effectively decou- 
ple by means of a powerlaw. Finally, the generic UV 
incompleteness of axion-electrodynamics puts an upper 
bound on the axion-photon interaction for small axion 
masses much in the same way as the Higgs boson mass 
is bounded from above in the standard model. 

Our article is organized as follows: in Sect. [TTJ we intro- 
duce axion electrodynamics in Minkowski and Euclidean 
space. Section [TTT| is devoted to the quantization of axion 
electrodynamics using the Wetterich equation, i.e., the 
RG flow equation for the 1PI generating functional. The 
RG flow is investigated in the theory space spanned by 
the axion mass and axion-photon coupling in Sect. |IV[ 
where the general structure of the flow is discussed, as 
well as explicit solutions are worked out. Phenomenolog- 
ical implications in the context of the QCD axion as well 
as for more general axion-like particles (ALPs) are sum- 
marized in Sect. |Vj We conclude in Sect. I VII and defer 
some technical details (App.lAl) as well as further con- 
ceptual considerations (App. [B]) to the appendices. 



II. AXION ELECTRODYNAMICS 

We consider an axion-photon field theory with clas- 
sical Lagrangian in Minkowski space (with metric rj — 
(+>---)) 

C = -l F ^F^-^d,ad^a-^m 2 a 2 -^gaF, u F^, (1) 

where m denotes the axion mass and g the axion-photon 
coupling. The latter has an inverse mass dimension. This 
Lagrangian - considered as the definition of an effective 
classical field theory - serves as a starting point for a 
variety of axion-photon phenomena such as the Primakoff 
effect [TT] or axion-photon oscillations giving rise to light- 
shining-through- wall signatures [16], polarimetry effects 
[IHIEO] or higher-frequency generation [3TJ. For the QCD 
axion, the mass and coupling parameters are related, g ~ 
m. Particular implications for the example of a QCD 
axion will be discussed in Sect. |Vj In the following, we 
take a different viewpoint and consider this Lagrangian 
as a starting point for an (effective) quantum field theory. 

Let us carefully perform the rotation to Euclidean 
space where the fluctuation calculation will ultimately 



be performed. We use the standard conventions 

x"\m = -ixi\B, A°\ M = -iA 4 \ E , <9°| M = «<9 4 |e, 
£|m = -£|e, (2) 

where |m/e marks the quantities in Minkowski/Euclidean 
space, respectively. As a consequence, we get, for 
instance, 9 M a9^a| M = — d^ad^a^, and F^ U F^\ M — 
FfiuFftul-E, as usual. The axion-photon coupling con- 
tains the structure aF^F^ = ^ae^^F^ F kX , which 
includes d°A 1 d 2 A 3 \M and index permutations thereof as 
building blocks. As only one time-like component ap- 
pears, the Minkowskian and Euclidean versions of the 
axion-photon couplings differ by a factor of i. A possible 
further but irrelevant factor of (—1) may or may not arise 
depending on the conventions for the Levi-Civita symbol. 
The resulting Euclidean Lagrangian finally reads 

£|e = ^F^F^ + ^d^ad^a+^m 2 a 2 + ^igaF^F^, (3) 

For a computation of photon fluctuations in the contin- 
uum, gauge fixing is necessary. As there are no mini- 
mally coupled charges in pure axion electrodynamics, we 
can perform a complete Coulomb- Weyl gauge-fixing pro- 
cedure, imposing the gauge conditions 

V • A = 0, A = 0, (4) 

such that only the two transversal degrees of freedom of 
the photon field remain. The corresponding gauge-fixing 
action in Minkowski space reads 

£gf = ^-(rT - n»n v )(drA v )(d"A x )(ri lt x - n K n x ) 
la 

+^(n^) 2 , (5) 

where = (1,0,0,0). The translation to Euclidean 
space using Eq. ([2| is straightforward. The two gauge- 
parameters a, j3 can be chosen independently. In the 
present work, we will exclusively consider a Landau- 
gauge limit, a,j3 — > 0, corresponding to an exact im- 
plementation of the gauge-fixing conditions. 

III. QUANTIZATION OF AXION 
ELECTRODYNAMICS 

As the axion-photon coupling has an inverse mass di- 
mension, i.e., [g] — —1, axion electrodynamics does not 
belong to the class of perturbatively renormalizable theo- 
ries. Without any further prerequisites, quantum fluctu- 
ations in this theory can still be dealt with consistently 
in the framework of effective field theories. This may 
potentially require the fixing of further physical parame- 
ters corresponding to higher-order operators. If the the- 
ory featured a non-Gaufiian UV fixed point with suitable 
properties, axion electrodynamics could even be asymp- 
totically safe [35] |33J and thus non-perturbatively renor- 
malizable. 
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A suitable framework to deal with quantum fluctua- 
tions in either scenario is the functional renormalization 
group. It allows us to study the renormalization flow 
of rather general classes of theories specified in terms of 
their degrees of freedom and their symmetries. For this, 
a convenient tool is the Wetterich equation [34] 



d t T k = ^Tr 



d,R k 



.(2) 



Rk 



d t = k 



dk' 



(6) 



for the effective average action which interpolates be- 
tween a microscopic or bare UV action 5a = and 
the full quantum effective action T = T^q. The effec- 
tive average action governs the dynamics of the field 
expectation values after having integrated out quantum 
fluctuations from a UV scale A down to the infrared scale 
k. The infrared regulator function Rk specifies the de- 
tails of the regularization of quantum fluctuations to be 
integrated out near an infrared momentum shell with mo- 

(2) 

mentum k. T k denotes the second functional derivative 
of the effective average action with respect to the fields, 
and the trace contains a summation/integration over all 
discrete/continuous indices, reducing to a momentum in- 
tegral in the simplest case. Thus the Wetterich equation 
is a one-loop equation from a technical point of view, 
but nevertheless includes effects at higher loop order in 
perturbation theory since it is the full non-perturbative 
propagator that enters the loop diagram. Since its deriva- 
tion does not rely on the existence of a small parameter, 
it is applicable also in the non-perturbative regime. For 



reviews of the functional RG see, e.g., [55] . 

In the present work, we study the renormalization flow 
of axion electrodynamics as parameterized by a class of 
(Euclidean) action functionals of the form 



d A x 



4 



z, 



Z F / d 4 x£ GF ,(7) 



a(x)F ilv {x)F llu {x) 



where we allow the mass and the coupling constant to be 
scale dependent. Also, the wave function renormaliza- 
tions Zpj a are implicitly assumed to be scale dependent. 
In principle, also the gauge-fixing parameters a, /? could 
run with k. However, we use the known fact that the 
Landau-gauge limit is a fixed point of the RG 36J . Fur- 
ther terms beyond this simple truncation of full axion- 
electrodynamics are induced by the operators present in 
our truncation, even if set to zero at a UV scale A, and 
are generically expected to couple back into the flow of 
the couplings considered here. 

The flow equation involves the second functional 

(2) 

derivative of the action T k , corresponding to the inverse 
propagator, which is matrix-valued in field space, 



r (2) _ 

1 k — 



(8) 



In momentum space, the corresponding components read 



T^aA 

1 A — 



5a(p)Sa(— q) 

S 2 T k 
8a{p)5A x {-q) 

5 2 r k 

5A K (p)Sa(-q) 



(2n) i S^(q-p)[m 2 +Z a q 2 }, 
= igke^ 15 5 vX A s {q - p)q l _ l (q 1 - _p 7 ) , 



*9k£ 



^ S S^A s (q-p)p^ Pl -q 7 ). 



(9a) 
(9b) 
(9c) 



1 kA 



6A K (p)6A x (-q) 

i9ke^ S 6^Ss^ Pl a(q-p) + (2i:) A 5^{q - p)Z F 



(q 2 S K \ - q K q\) + (S KfJ . - n K n li ) < ^ ( ^{8 y x ~ n v ti\) + 



(9d) 



Note that T aa and the second term of T AA denote the field- independent inverse propagators, whereas T aA , T Aa and 
the first term of T AA correspond to vertices, since they contain powers of the external fields. The regulator R k is 
chosen diagonal in field space with the components, 



RT{q) = Z a q 2 r(q 2 /k 2 ), R AA ^(q) = Z F q 



S K \ - 



aq z p q l 



The photon propagator in Landau gauge then takes the form 

1 



r{q 2 lk 2 \ 
(10) 



(TW+R k )-') AA (q)\ a = = 



1 / 2 \ 71 ' 1 l I 

o K \ ^ ; ^ {q K ,Q\ + q n K n\ ) + — ; ^ (n R q\ + n x q K ) , 



W (1 + rjcf A 2 )) \ kX q 2 -(n- g) 2 ^ + Q """^ + q 2 - (n ■ qf 



(11) 
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r 



where r(y) denotes a regulator shape function, specify- 
ing the details of the regularization scheme, see below. 
Evaluating the Wetterich equation and determining the 
flow of Tfe from a UV cutoff A down to k = corresponds 
to a quantization of axion electrodynamics in the path- 
integral framework. If a cutoff A — > 00 limit existed, ax- 
ion electrodynamics would even be UV complete. In the 
following, we determine the flow of Tk truncated down to 
the form of Eq. ([7]) as a nonperturbative approximation 
to the full effective action. 



pseudo-Goldstone boson of a spontaneously broken sym- 
metry. In other words, an axion mass term or a higher- 
order self-interaction potential cannot be generated di- 
rectly from the axion-photon vertex. A simple proof of 
Eq. (16) can be based on the observation that both the 



axion mass term as well as the axion-photon coupling 
are operators which are nonzero for constant axion fields, 
a = a c . Hence, a projection of the right-hand side of the 
flow equation onto a = a c suffices in order to extract the 
flow of these operators. However, for a constant axion, 
the axion-photon coupling can be written as 



IV. RG FLOW OF AXION 
ELECTRODYNAMICS 

We are interested in the RG flow of all scale-dependent 
parameters of Tk in Eq. Let us start with the flow 
of the wave function renormalizations Z a and Zp . Their 
flow can be extracted by projecting the right-hand side 
of the Wetterich equation onto the corresponding kinetic 
operators: 



d t Z a = 
d t Z F = 



dq 2 
4 d 



d 4 p 



s 2 d t r k 

5a(p)5a(—q) 



d pn K n x 



a,A,q-^0 

d t T k 



5A K (p)SA x (-q) 



(12) 



(13) 

a,A,q— >0 



where f2 denotes the spacetime volume. Similar projec- 
tion prescriptions also exist for the flow of the (unrenor- 
malized or bare) axion mass and axion-photon coupling, 



d t m k 
d t gk 



1 

ft 
1 i 



d 4 p 



Vl 24 a " dr a dqp 



H*- k 



8a(p)Sa(—q) 
d d 



(14) 



a,A,q— ¥0 



d 4 P 



HL k 



Sa(p)SA K (r)6A x (-q) 



.(15) 



a,A,q— >0 



It is straightforward to verify that - within our truncation 
the RG flow of the latter quantities vanishes exactly, 



d t ml 



0, d t g k = 0. 



(16) 



The reason for this nonrenormalization property is obvi- 
ous from the structure of the axion-photon vertex. From 



the first term of Eq. ( 9d I , we observe that the vertex with 



an external axion field has a nontrivial momentum struc- 
ture ~ qpa(q — p). Together with the anti-symmetric 
Lorentz structure this vertex can only generate opera- 
tors containing derivatives of the axion field ~ da. This 
also ensures that no non-derivative axion-fermion cou- 
plings can be generated by integrating out photon fluc- 
tuations, as is in accordance with the axion being the 



d X Oj c Fui/Fny 2ft c 



d 4 zd M (a^ 



(17) 



i.e., as a total derivative of an abelian Chern-Simons cur- 
rent. Therefore, it can be eliminated from the action and 
thus cannot contribute to the renormalization flow of the 
nonderivative axion terms of the theoryj^] 

The nonrenormalization property extends to a full ax- 
ion potential 



d t V k (a 2 ) = 0, 



(18) 



which holds as long as the potential is purely mass-like 
at some initial scale A, V^(a 2 ) = \m\a 2 . Of course, 
as soon as axion self-interactions are present at some 
scale, all parameters of the potential get renormalized 
through these self-interactions. In other words, the pure 
axion mass term and the axion-photon coupling are par- 
tial fixed points of the RG flow of axion electrodynamics 
in terms of bare quantities. 

The only renormalization effects in our truncation 
therefore arise from the running wave function renormal- 
izations. The physical observables can be expressed in 
terms of the (dimensionful) renormalized couplings, given 
by 



si 



ZpZ a 



(19) 



In order to investigate the structure of the RG flow and 
in particular to search for fixed points at which the the- 
ory becomes scale-free, it is convenient to introduce the 
renormalized dimensionless axion mass and axion-photon 
coupling, 



k 2 Z a 



k 2 ' 



.9 2 = 



fkk 2 
Z\Z a 



2 ; 2 
= 5R fc 



(20) 



A similar observation was made in |37| in order to inversely ar- 
gue that the renormalization flow of the topological charge in 
a non-abelian gauge theory can be properly formulated if the 
topological charge 8 is temporarily considered as a spacetime 
dependent field. 
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a 



(b) 




Figure 1. Schematic representation of the diagrammatic con- 
tributions to r]p (a) and r\ a (b). Dashed lines denote axions, 
curly lines photons. As implied by the exact flow equation, 
all internal propagators and vertices are considered as fully 
dressed at the scale k. Appropriate insertions of dtRk at one 
of the internal lines in each diagram and a corresponding sum 
over insertions is understood implicitly. The contribution to 
r/F is thus to be understood as two distinct diagrams, associ- 
ated with the contributions ~ rj a and ~ t)f in Eq. (251. 



Defining the anomalous dimensions of the axion and the 
photon field 



rj a = -d t In Z a , r/F = -d t In Z F , 



(21) 



the ft functions of mass and coupling can be written as 
d t g 2 =P g 2 = (2 + 2r ?F + ry Q ) ff 2 , (22) 



d t m 2 = /3 m 2 = (rj a - 2)m 2 . 



(23) 



From the projections rtl2|) and (13 



of the flow onto the 
kinetic operators, the anomalous dimensions can be ex- 
tracted: 



la 



Vf 



6(4tt) 2 

„2 



(2-¥) , (2 



Vf ■ 
4 . 



6(4tt) 2 V(l + m 2 ) 2 



(24) 



(25) 



For these specific forms, we have used the linear regu- 
lator shape function r(y) = (i — 1)9(1 — y) [38 . The 
corresponding results for arbitrary shape functions are 
given in Appendix [3] These are the central results of the 
present paper, diagrammatically represented in Fig. [T] 

As axion electrodynamics belongs to the perturba- 
tively non-renormalizable theories according to naive 
power-counting, the above-given ft or 77 functions do not 
exhibit the same degree of universality as ft functions 
of marginal couplings in renormalizable theories. This 
implies that typically even the leading-order (i.e., "one- 
loop") ft function coefficients are scheme dependent. In 
the present case, we can change the numerical value of 



the prefactors in Eqs. (24) and (25) by varying the reg- 
ulator shape function rjy). Nevertheless, the sign of the 
prefactors cannot be changed for admissible regulators as 
is visible from the explicit representations given in Ap- 
pendix [X] Furthermore, even for theories exhibiting this 
type of non-unversality in their ft functions, the existence 
of fixed points and the critical exponents determining the 
universality class of the fixed point are universal. 

Let us first discuss the resulting ft functions in various 
simple limits. At weak coupling, the anomalous dimen- 
sions are small, rj ~ g 2 . This implies that the anomalous 



dimensions on the right-hand sides of Eqs. (24) and (25), 
signaling a typical RG "improvement" , can be neglected. 
Further assuming a heavy-axion limit, m 3> 1, corre- 
sponding to a renormalized axion mass being larger than 
a given scale under consideration, m 2 ^ 3> k 2 , the axions 
decouple, yielding 

Q 2 

rip^O, ria^^r-^, form>l,.g<l. (26) 
3(47r)^ 

The axion anomalous dimension remains finite due to 
the fact that the massless photons never strictly decouple 
from the flow. In the present weak-coupling heavy-axion 
limit, the remaining flow of the coupling reduces to 



d t g 2 = 2 



9 



3(4tt) : 



9 



for m> l,g < 1, (27) 



which can be straightforwardly integrated from a high 
UV scale A to some low scale k. In terms of the dimen- 
sionful renormalized coupling, the solution reads 



gUk) 



ffl(A) 



1 



6(4tt) 



(A2-P) 5 |(A)' 



for to 3> 1 , g <C 1 , 
(28) 

where <?|,(A) denotes the initial coupling value at the UV 
cutoff. This value has to satisfy the weak coupling condi- 
tion g 2 i (A)A 2 -C 1 • In the present weak-coupling heavy- 
axion limit, we conclude that the photon-axion coupling 
undergoes a finite renormalization even in the deep in- 
frared (IR) limit, despite the fact that photonic fluctua- 
tions never strictly decouple in the deep IR. The photons 
still effectively decouple, as their low-momentum contri- 
butions to the flow vanish according to the powerlaw 
~ A: 2 <7r,(A) for k — > 0. The same conclusion holds for 
the axion mass. Inserting the coupling solution (28) into 
the flow for the mass ( 23 1 , the solution for the renormal- 
ized mass reads 



i(A) 



l + 6^F(A 2 -fc 2 )ff|(A)' 



for m»l,j«l, 
(29) 

Axion electrodynamics therefore exhibits a remarkable 
IR stability. Quantitatively, both the coupling as well as 
the axion mass run to slightly smaller values towards the 
infrared in the present limit. 

Let us now turn to the massless axion limit. As 
is obvious from Eq. (23), the massless theory is an 
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RG fixed point, m* = 0. This follows from the non- 
renormalization of the axion potential, which also im- 
plies that the axion will presumably not exhibit the fine- 
tuning problem typically associated with massive scalars. 
As long as r\ a < 2, this fixed-point is IR attractive. For 
weak coupling, the anomalous dimensions reduce to 



2 g 2 g 2 
VF ^3{4^' Va ~*3(4^' form = '9« 1 ' ( 3 °) 



Integrating the flow analogously to Eq. ( 28 1 leads us to 



yielding the coupling flow 



<hu L = ■■[ 2 + 5^|^ J g 2 , for m = 0, g « 1. (31) 



5|(A) 



1 



6(4tt) 



r (A2-fc2) 5 |( A) ' 



for m — 0, g <C 1. 

(32) 

The conclusion is similar to the heavy-axion case: even 
though there is no decoupling of any massive modes, ax- 
ion electrodynamics shows a remarkable IR stability. The 
fluctuations of the massless degrees of freedom induce 
only a finite renormalization of the axion-photon cou- 
pling yielding smaller couplings towards the IR. 

In the intermediate region for finite but not too heavy 
masses, m* = remains an IR attractive fixed point 
for weak coupling, but this massless point is so weakly 
attractive that decoupling of the axions typically sets in 
first and the flow ends up in the heavy-axion limit. 

Let us now turn to arbitrary values of the coupling. 
For this purpose, we solve Eqs. (24) and (25) for r) a ,rjF 
and insert the result into Eqs. (22) and (23), yielding the 



obviously nonperturbative flow equations 



d t g 2 = = 2g : 



13g 4 - 384ttV(21 + 17m 2 + 4m 4 ) - 147456vr 4 (l + m 2 ) 2 



g 4 - 384ttV(1 



147456tt 4 (1 



d t m 2 = (3 rl 



6m 



1287rV(4m 4 + 7m 2 + 3) - 49152tt 4 (1 



,2\2 



384ttV(1 



147456tt 4 (1 



,2^2 



(33) 
(34) 



In the heavy-axion limit the decoupling of the axion fluc- 
tuations still simplifies the system considerably, since the 
suppression of rip{m — > oo) — > still persists, leading 
again to 



d t g 2 = 2 



9' 
3(4tt) 2 



for m 3> 1, 



(35) 



as in Eq. (27). Also the corresponding mass flow is equiv- 



alent to that of the heavy-axion weak-coupling limit. The 
resulting flow exhibits no signature of a fixed point apart 
from the Gaufiian one at g — 0. On the contrary, the /3 
function for the running coupling is somewhat similar to 
that of one-loop QED. The integrated flow will therefore 
give rise to a Landau pole of the running coupling at high 
scales. Fixing the renormalized dimensionful coupling to 
some value at the scale k = 0, i.e., <?r(/c = 0) = <?p , the 
scale of the Landau pole /xl, where g^iph) oo, can 
directly be read off from Eq. |28| , 



ML 



9ro 



31 

9ro' 



(36) 



which is an order of magnitude larger than the inverse 
coupling in the deep IR. For this estimate to hold, the 
dimensionless mass has to be large during the whole flow. 
With regard to Eq. (29), this is true as long as 1 < 



m 2 = rn/^/k 2 is satisfied for all k. This is always true if 



m|(fc 



1 (this criterion could even be relaxed 



a bit). Going towards larger scales, mp{(A:) increases with 
k and diverges even at k = fi^. 

A different structure becomes visible in the massless 
limit for arbitrary coupling. Whereas m* = still is a 
fixed point of the RG, the nonperturbative flow of the 
coupling becomes 



d t g 2 = 2g 



2 13 5 4 - 8064ttV - 1474567T 4 
g A - 384ttV - 1474567T 4 



for m = 0. 



(37) 

Starting from small coupling values, the (3 g 2 function runs 
into a singularity (/3 g 2 — > oo) for 



<? s 2 ; = 64tt 2 (3 + 3 V5) , for m = 0, 



(38) 



where the denominator changes sign. For even larger 
couplings, the /3 g 2 function returns from — oo and has a 
zero at 



13 



for m = 0. 



(39) 



For g 2 — y oo, the /3 g 2 function approaches the simple form 
d t g 2 ~ 26g 2 , for g 2 -> oo, m = 0. (40) 
The massless f3 g 2 function is plotted in Fig. [2] 
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Figure 2. The /3 function /3 g 2 in the massless limit clearly 
exhibits an IR-attractive GFP, a singularity at g = fusing, and 
an asymptotic behavior ~ 26g 2 . 



Fixing the renormalized dimcnsionful coupling in the 
deep IR as above, <? R (0) = 5roj we can integrate the flow 
until the running coupling hits the singularity of the /3 g 2 
function. The corresponding scale fi'^ can be viewed as 
an "RG-improved" Landau pole, giving a nonperturba- 
tive estimate of the scale of maximum UV extension of 
quantum axion electrodynamics. From a numerical inte- 
gration, we obtain 



14.844 

ffRO 



for m = 0. 



(41) 



This is of the same order of magnitude as the scale 
found above for the heavy-mass limit. As a rough esti- 
mate of the regularization-scheme dependence of this re- 
sult, let us note that in the case of an exponential shape 
function r{y) = (exp(y) - l)" 1 we obtain // Lexp ~ 

For finite values of the axion mass, the flow of the 
coupling generically exhibits the same features as for the 
massless case. The position of the singularity in j3 g 2 is 
shifted towards larger values of the coupling for increas- 
ing mass. It is straightforward to verify that the j3 func- 
tions ( 33 1 , ( 34 ) do not support a nontrivial fixed point 
for physically admissible positive values of g 2 and m 2 . At 
the same time, the mass flow exhibits the same singular- 
ity in /3 m 2 as the denominators in Eqs. (33) and (34) are 
identical. 

We conclude that generic flows in the physically admis- 
sible parameter space cannot be extended beyond a scale 
of maximum UV extension. Phenomenological implica- 
tions of the existence of such a maximum UV scale will 
be discussed in the next section. Whether or not sub- 
stantial extensions of our truncation are able to modify 
this conclusion is hard to predict. For instance, an inclu- 
sion of axion self-interactions would add another sector 
to the theory which generically suffers from a triviality 
problem and is thus not expected to change our conclu- 
sions. Qualitative modifications of our j3 functions could 
potentially arise from momentum-dependent axion self- 
interactions. These are induced by photon fluctuations 
even within our truncation, and can couple into the flow 



of r] a , thus also contributing to the j3 functions for the 
axion mass and axion-photon coupling, as well as further 
momentum-dependent axion self interactions. 

For the remainder of this section, let us concentrate 
on an oddity of the present flow. There is in fact one 
exceptional RG trajectory which can be extended to all 
scales. This trajectory requires the singularity induced 
by the denominator of the (3 function to be canceled by a 
zero of the numerator. Remarkably, there exists a mass 
and coupling value in the physically admissible region, 
where the singularities in both the mass and the coupling 
flow are canceled. This exceptional point in theory space 
is given by 



= ~(V5-1), 



Q 2 

CXC 

6(4tt)< 



= 2(3 



(42) 



The RG trajectory which passes through this point, say 
at a scale A exc , has a standard IR behavior exhibiting 
the typical decrease towards smaller renormalized mass 
and couplings. Towards the UV, the f3 functions do not 
exhibit a fixed point but approach the simple form 



d t g 2 ~ 26g 2 , d t m 2 



6m 



(43) 



which resembles a pure dimensional scaling with large 
anomalous dimensions. In fact, this behavior corresponds 
to rj a — > 8, and i]p — > 8. It implies that both di- 
mensionless and dimcnsionful renormalized couplings in- 
crease strongly towards the UV without hitting a Landau 
pole singularity. Instead the growth of the couplings re- 
mains controlled on all scales and approaches infinity at 
infinite UV cutoff scale. 

As expected, the scaling dimensions are non-universal, 
and in fact show a considerable regulator dependence: 
Employing an exponential shape function yields dtg 2 — 
^g 2 and d t m 2 ~ fm 2 , corresponding to r] a — 4.5 = r/p. 

An important consequence of this exceptional flow is 
that the physical IR values are completely fixed in terms 
of the scale A oxc . Numerically integrating the flow to- 
wards k — > yields for the linear regulator 



g R (k = 0) = g m 



20.36 

AqxC 

m R (fc = 0) = toro ^ 0.4657A e 



(44) 



(45) 



implying the dimensionless combination <?roTOro — 
9.484. Of course, our observation of this exceptional tra- 
jectory requires a critical discussion: The fact that the 
anomalous dimensions become comparatively large may 
be interpreted as a signature that explicit momentum- 
dependencies of the propatators and vertices become im- 
portant at larger coupling. If so, the exceptional tra- 
jectory might just be an artefact of our truncation which 
assumes tree-level-type propagators and vertices. On the 
other hand, we cannot exclude the possibility that this 
exceptional trajectory is a simple projection of a legiti- 
mate UV-extendable trajectory in the full theory space. 



If so, the observed exceptionality could reflect the restric- 
tions on the physical parameters induced by the true UV 
behavior (potentially controlled by a UV fixed point). 
Our truncated RG flow could then be quantitatively re- 
liable below A oxc . The strong regulator dependence ob- 
served above should be read as a hint that if such an ex- 
ceptional trajectory exists within the full theory space, 
it might still change considerably when operators beyond 
our present truncation are taken into account. In partic- 
ular, the large anomalous dimensions seem to call for the 
inclusion of higher-derivative operators. 

In the rather speculative case, that the trajectory ex- 
ists with the same qualitative features in the full the- 
ory, we would have discovered a (to our knowledge) first 
example of a UV complete theory with a high energy 
behavior that is controlled by neither a fixed point nor 
a limit cycle, see, e.g., [39 . Though this theory would 
not fall into the class of asymptotically safe systems due 
to the lack of a UV fixed point, it would be asymptot- 
ically controllable. The number of physical parameters 
of such systems would then correspond to the dimen- 
sionality of the exceptional manifold, i.e., the analog of 



Eq. (42 1 including all possible further couplings. We em- 



phasize that in this speculative case also further proper- 
ties required for a legitimate field theory such as unitar- 
ity would have to be critically examined. In particular, 
a large positive anomalous dimension, as observed here, 
implies a strongly UV-divergent propagator ~ (p 2 ) 1- '^ 2 , 
which might result in cross-sections increasing as a large 
power of the momentum. Whether such behavior can be 
reconciled with requirements such as perturbative uni- 
tarity within standard quantum field theory remains to 
be investigated. 



V. PHENOMENOLOGICAL IMPLICATIONS 

A. QCD axion 

Axion electrodynamics occurs naturally as a low- 
energy effective theory in the context of the Peccei-Quinn 
solution of the strong CP problem. The QCD axion de- 
velops a generic two-photon coupling, induced by its mix- 
ing with the mesonic 7r°, 77 and 7/ degrees of freedom. 
Effective axion electrodynamics therefore arises during 
the chiral phase transition of QCD at around a typical 
QCD scale which we choose to be Aqcd — 1 GeV. At 
that scale, the relation between axion coupling and mass 
is essentially fixed by the corresponding pion scales, 



Sr(Aqcd) = C 



™r(Aqcd) 

f-K 



(46) 



where and f n denote the pion mass and decay con- 
stant respectively, and C is a numerical dimensionless 
constant that depends of the microscopic axion model 
(typically defined in the context of a grand unified sce- 
nario). For instance, for the KSVZ axion [3J, the cou- 
pling mass relation yields £(r(Aqcd) — ^2"1r(Aqcd)- 



Other axion models typically lie within an order of mag- 
nitude near this relation [4J [40] . 

Axion searches in experiments or astrophysi- 
cal/cosmological observations do actually not test 
the parameters occurring in Eq. ( 46 ) directly, as they do 



not operate near the QCD scale, but typically at much 
lower scales fc bs- They range from ^keV momentum 
scales of stellar evolution theory of horizontal-branch 
stars, via ~eV scales for solar energy loss or direct 
helioscope observations, to ~ fieV scales in light-shining- 
through walls experiments. In other words, a proper 
comparison of such observational results with Eq. ( 46 ) 



requires to take the finite renormalization of axion elec- 
trodynamics between Aqcd and fc bs m t° account. For 
simplicity, we estimate the maximum renormalization 
effect by choosing fc b s = 0. In order to illustrate our 
findings, we plot in Fig. [3] the renormalized ratio 



cr 



gR(fcobs)/wR(fc b s ) 

gn (Aqcd ) /m R ( A QC d ) 



(47) 



inspired by Eq. ( 46 ) as well as the ratio of the product 



PR 



g R (fcobs)ra R (fc obs ) 
Sr(Aqcd)™.r(Aqcd) 



(48) 



as dark-blue (dashed) and purple (solid) lines, respec- 
tively. The ratios are plotted as a function of the loga- 
rithm of the initial coupling 5r(Aqcd) in units of GeV. 
For small couplings, the renormalization of both coupling 
and mass remain unobservably small implying that both 
ratios cr,pr — > 1. Only for larger couplings, the renor- 
malization of both towards smaller values becomes visi- 
ble, implying a significant decrease of pr. Most impor- 
tantly, the proportionality between mass and coupling 
remains essentially unaffected. Numerically, the ratio cr 
changes only on the 10 -5 level. 

We conclude that exclusion bounds derived from the 
non-observation of axion effects are not modified by 
renormalization effects of axion electrodynamics along 
the lines of constant <7r/Wr. As the physically rel- 
evant parameter space for the QCD axion lies below 
<7r < 10 _9 (GeV) _1 , we conclude from Fig. [3] that the 
renormalization effects within the effective theory of ax- 
ion electrodynamics are completely irrelevant. 

The existence of a scale of maximum UV extension 
/il also seems to be of no relevance for the QCD ax- 
ion: for instance, for an axion coupling in the range of 
g ~ O(10 _16 )GeV where the axion could also provide a 
substantial part of the dark matter of the universe, the 
scale /jl would lie well beyond a typical GUT scale. In 
turn, this demonstrates that axion electrodynamics is a 
consistent effective (quantum) field theory at low energies 
within the QCD axion scenario. 

For completeness, let us mention that the QCD axion 
would be on the exceptional trajectory found above, if 
Sr(Aqcd) ^ 1.9 Ge^ 1 andm R (A QCD ) - 5.0 GeV. This 
is, of course, beyond the parameter regime conventionally 
considered for the QCD axion. 
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Figure 3. Plot of the renormalized coupling ratios cr, of 
(47 1 (blue/dashed line) and pn of (481 (purple/solid line) as 



a function of the logarithm of the initial axion-photon cou- 
pling (/r(Aqcd) in units of GeV. For small couplings, the 
renormalization of both coupling and mass remain unobserv- 
ably small. Only for larger couplings, the renormalization 
of both towards smaller values becomes visible (purple/solid 
line), while the proportionality between mass and coupling 
remains essentially unaffected (blue/dashed line). 



B. General axion-like particles (ALPs) 

Apart from the QCD axion, massive scalar bosons 
may arise as pseudo-Nambu-Golstone bosons in a vari- 
ety of contexts [3 [8]- In case of a pseudo-scalar bo- 
son, the emergence of a coupling to photons of axion- 
electrodynamics-type is a natural consequence. Since the 
coupling-mass relation in this general case is not neces- 
sarily fixed as in Eq. ( 46 1 , a more general parameter space 



for such axion-like-particles can be investigated. 

A priori, the mass-coupling space seems rather unre- 
stricted. Depending on the mass-generation mechanism, 
the ALP mass can vary strongly from very light masses 
arising from anomalous breakings as in the QCD-axion 
case, to explicit masses of the order of the high scale Atjv 
(such as the Fermi, GUT, or Planck scale). Depending 
on the value of the coupling, the low-energy mass may be 
very different from the high-energy mass due to renormal- 



izations of the type ( 29 ) , but, in general, the renormalized 
low-energy mass still remains rather unconstrained. 

This can be different for the possible values of the ef- 
fective low-energy ALP-photon coupling. If the effective 
low-energy theory still exhibits a scale of maximum UV 
extension /j,^ as in Eqs. (36) or (41 ), then that scale must 



necessarily be higher than the microscopic high energy 
scale Atjv- Otherwise, the effective low-energy theory 
could not possibly arise from that unknown microscopic 
theory. As a consequence, the low-energy coupling is 
bounded from above due to the renormalization flow in 
the effective theory, 



maxsr R (fc -> 0) = gno,r, 



Nr < N K 
Ml 



A 



uv 



(49) 



where Nr is a number that depends on the ALP mass 
as well as on further degrees of freedom in the effective 



low-energy theory. In the case that this effective theory 
is well-approximated by axion electrodynamics, we found 
that 2Vr = O(10) (A/r ~ 31 for large masses and Ar, — 
14.8 for the massless case). 

Our conclusion, that g^o should be suppressed by the 
high scale Atjv looks rather trivial, as it seems to fol- 
low standard power-counting arguments for a higher- 
dimensional operator. However, we stress that the state- 
ment is actually stronger: even for unusually enhanced 
couplings at the high scale (invalidating naive power- 
counting), axion electrodynamics as an effective theory 
ensures that the low-energy coupling is supressed by 
renormalization effects and obeys the bound ( 49 ) . 



We expect the precise number JVr. to be modified by 
further degrees of freedom, such as the standard-model 
fermions contributing to the flow at higher scales. As 



long as they leave the anomalous dimension r) a in Eq. ( 24 ) 



positive, we expect a mere quantitative influence on Nr. 
Of course, these conclusions no longer hold, if the ad- 
ditional degrees of freedom render the effective theory 
asymptotically free or safe, or if the system sits on the 
exceptional trajectory as discussed above. 



VI. CONCLUSIONS 

We have analyzed the renormalization flow of axion 
electrodynamics considered as an effective quantum field 
theory. From a field-theory viewpoint, we have revealed 
several interesting properties: non-renormalization prop- 
erties protect the flow of the axion mass and axion- 
photon coupling, such that within our truncation the 
flow of the corresponding renormalized quantities is solely 
determined by the axion and photon anomalous dimen- 
sions. These non-renormalization properties are in line 
with but go beyond the fact that the axion can be under- 
stood as a pseudo-Nambu-Goldstone boson of a broken 
Peccei-Quinn symmetry. 

Towards the infrared, the flow remains well controlled 
even in the presence of massless degrees of freedom. Even 
though massless photons strictly speaking never decou- 
ple, their contribution to the RG flow effectively decou- 
ples leading to finite and predictable IR observables. By 
contrast, the UV of axion electrodynamics exhibits a triv- 
iality problem somewhat similar to 4 theory: insisting 
on sending the cutoff A — > oo is only possible for the 
free theory g — > 0. From a more physical viewpoint, fix- 
ing the coupling to a finite value at a finite scale implies 
that axion electrodynamics can be treated as a quantum 
field theory only up to a scale of maximum UV extension 
/.tL- The value of this scale is, of course, not universal. 
For our regularization scheme, this scale is about an or- 
der of magnitude larger than the inverse IR coupling and 
depends weakly on the axion mass. 

A behavior different from this generic case is only 
found on an exceptional RG trajectory which remains 
free of singularities on all finite scales. From within our 
truncated RG flow, it is difficult to decide whether this 
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trajectory is a mere artefact of the truncated theory space 
or a remnant of a valid trajectory of an interacting UV- 
controlled theory. In the latter case, it could be the first 
example of a predictive and consistent theory without 
being associated with an obvious UV fixed point. At the 
present level of approximation, however, we consider this 
exceptional trajectory as an oddity, the status, physical 
relevance and consistency of which still has to be care- 
fully examined. 

For the QCD axion, our findings demonstrate that the 
typical proportionality between axion-photon coupling 
and axion mass is not relevantly renormalized by low- 
energy axion-photon fluctuations. In fact, possible renor- 
malizations of mass and coupling largely cancel out of the 
proportionality relation even at stronger coupling. As a 
consequence, phenomenological bounds on these param- 
eters remain essentially unaffected by the RG flow. By 
contrast, the absolute values of mass or coupling can un- 
dergo a sizeable renormalization, however, the required 
coupling strength is not part of the natural QCD axion 
regime. 

For more general axion-like particles the existence of 
a generic maximum scale of UV extension induces an 
upper bound on possible values of the axion-photon cou- 
pling. Standard lines of perturbative reasoning |44j sug- 
gest that the renormalized coupling should be of order 
5r ~ 1/Atjvj where Auv denotes the microscopic scale 
where the axion sector is coupled to the standard model 
particles. Our RG study now demonstrates that cou- 
plings of that size are not only natural, but are in fact 
bounded by <?r < O(10)/Auv due to renormalization ef- 
fects in the axion-photon sector. In view of the rather 



unconstrained ALP parameter space at large masses and 
comparatively large couplings (see, e.g., the compilation 
in [2]), our bound could become of relevance for ALP 
searches at hadron colliders above tor > lGeV and cou- 
plings above gn > 10 _3 /(GeV) -1 . 
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Appendix A: RG flows for general regulator shape 
functions 



For generality, we present the flows of the wave func- 
tion renormalizations for general regulator shape func- 
tions r(y). For the flow of the axion wave function, we 
obtain 



d t Z a 



2(2tt)' 3 



-VFV k (£) + d t r k (g) 

dp 



Z 2 F P 2 ( 



l + r fe 



The flow of the photon-wave function is given by 



(Al) 



d t Z F 
2 



d A p 



?aVaP 2 r k (£)+Z a p 2 d t r k (£) 



ZaP 2 



Zr 



-VFT k (£) +d t r k (jj) 



ZaP 2 



1 + 



(A2) 



Inserting the linear regulator shape function r(y) = 
(i — 1)6(1— y), the momentum integrals can be performed 
analytically. The results can be expressed in terms of the 



anomalous dimensions and are given by Eqs. (24) and 
pi 



Appendix B: Euclidean axion electrodynamics with 
imaginary coupling 

In the main text, we have emphasized that the phys- 
ically admissible parameter space is constrained to pos- 
itive masses and couplings (squared), g 2 ,m 2 > 0. For 
to 2 < and in absence of any further axion potential, 
the Euclidean action is unbounded from below along the 



direction of large axion field amplitude. Of course, this 
could be cured by adding a stabilizing potential, but this 
route will not be followed in this work. 



A more interesting case is provided by the case of 
imaginary axion-photon couplings, g 2 < 0. In this 
case, the Euclidean action would violate Osterwalder- 
Schrader reflection positivity, such that a corresponding 
Minkowskian theory can be expected to violate unitarity. 
Still, the Euclidean theory could be regarded as a valid 
field theory description of some suitable statistical sys- 
tem. In this case, the action would be stable along the 
axion- amplitude direction, but unstable towards the for- 
mation of large electromagnetic fields with large E ■ B. 
This instability could be cured by higher photon self- 
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interaction such as 

Cf = h{F^F^f + f 2 {F llv F llv f, (Bl) 

with positive constants /i^ carrying a mass dimension of 
—4. Actions of this type are familiar from fluctuation- 
induced nonlinear QED contributions to electrodynam- 
ics Also within axion electrodynamics, we ex- 



pect these contributions to be generated by mixed axion- 
photon fluctuations (some properties of such amplitudes 
have, for instance, been studied in [13]). In the following, 
we simply assume that these terms are suitably gener- 
ated either within axion electrodynamics or provided by 
an exterior sector coupling to photons. Then, also the 
parameter region where g 2 < can become physically 
admissible. 

In that case, the /3 function for g = —ig is given by 



2 13g 4 + 384tt 2 <? 2 (21 + 17m 2 + 4m 4 ) - 147456tt 4 (1 + m 2 ) 2 
9 5 4 + 384tt 2 5 2 (1 + m 2 ) - 147456tt 4 (1 + m 2 ) 2 ' ^ ' 



The mass /3 function now takes the form 



2 -g 4 - 1287r 2 5 2 (4m 4 + 7m 2 + 3) - 49152tt 4 (1 + m 2 ) 2 
m -g 4 - 384tt 2 5 2 (1 + m 2 ) + 147456tt 4 (1 + m 2 ) 2 ' ^ ' 



fit 



-150 




Figure 4. The /? function for the coupling g in the limit 
777 — ^ clearly shows a UV attractive non-Gaufiian fixed point 
at g — g* = 13.25, as well as an infrared attractive GauBian 
fixed point. 



In this system, we find a UV attractive fixed point at 
(g,m) ~ (13.25,0). The critical exponents, defined as 

eigenvalues of the stability matrix -g^- with — (<?,m), 
multiplied by an additional negative sign, are (^1,^2) = 
(2.16,1.20). Thus this fixed point is UV attractive in 
both directions. 

Clearly the system also admits a GauBian fixed point 
with critical exponents given by the canonical dimen- 
sions. Accordingly this fixed point is IR attractive in 
the coupling. 

We conclude that, setting the axion mass to zero, the 
system admits the construction of a complete RG tra- 
jectory, extending from g = at k — ¥ to g = for 
k — > 00, see fig. [4] Thus the axion-photon system with an 
imaginary coupling has the potential to provide a simple 
example of an asymptotically safe quantum field theory, 
albeit without an immediate physical application. 
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